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ABSTRACT 

Resonant relaxation (RR) is a rapid relaxation process that operates in the nearly-Keplerian potential near 
a massive black hole (MBH). RR dominates the dynamics of compact remnants that inspiral into a MBH and 
emit gravitational waves (extreme mass ratio inspiral events, EMRJs). RR can either increase the EMRI rate, 
or strongly suppress it, depending on its still poorly-determined efficiency. We use small-scale Newtonian TV- 
body simulations to measure the RR efficiency and to explore its possible dependence on the stellar number 
density profile around the MBH, and the mass-ratio between the MBH and a star (a single-mass stellar popu- 
lation is assumed). We develop an efficient and robust procedure for detecting and measuring RR in iV-body 
simulations. We present a suite of simulations with a range of stellar density profiles and mass-ratios, and mea- 
sure the mean RR efficiency in the near-Keplerian limit. We do not find a statistically significant dependence 
on the density profile or the mass-ratio. Our numerical determination of the RR efficiency in the Newtonian, 
single-mass population approximations, suggests that RR will likely enhance the EMRI rate by a factor of a 
few over the rates predicted assuming only slow stochastic two-body relaxation. 
Subject headings: black hole physics — galaxies: nuclei — stars: kinematics and dynamics — gravitational waves 



1. INTRODUCTION 

Dynamical relaxation processes near massive black holes 
(MBH) in galactic centers affect the rates of strong stellar in- 
teractions with the MBH, such as tidal disruption, tidal dissi- 
pation , or gravitational wave (GW) emission (e.g. lAlexanded 
l2005h . These relaxation processes may also be reflected by 
the d ynamical properties of the d ifferent stellar populations 
there dHop man & Alexandei*2006'), as observed in t he Galac- 
tic Center (Genzel et al. 2000; Pau mard et ani2006h . Of par- 
ticular importance, in anticipation of the planned Laser In- 
terferometer Space Antenna (LISA) GW detector, is to un- 
derstand the role of relaxation in regulating the rate of GW 
emission events from compact remnants undergoing quasi- 
periodic extreme mass ratio inspiral (EMRI) into MBHs. 

Two-body relaxation, or non-coherent relaxation (NR), is 
inherent to any discrete large-N system, due to the cumula- 
tive effect of uncorrected two-body encounters. These cause 
the orbital energy E and the angular momentum J to change 
in a random-walk fashion (ex \/t) on the typically long NR 
timescale Tnr. In contrast, when the gravitational poten- 
tial has approximate symmetries that restrict orbital evolu- 
tion (e.g. fixed ellipses in a Keplerian potential; fixed orbital 
planes in a spherical potential), the perturbations on a test 
star are no longer random, but correlated, leading to coherent 
(c>c t) torquing of J on short timescales, while the symmetries 
hold. Over longer times, this results in reso nant relaxation 
(RR) (iRauch & Tremainell996t[Rauch & Ing alls 1998; ' ^n\ . 
a rapid random walk of J on the typically short RR timescale 
T-RR <C Tnr. RR in a near-Keplerian potential can change 
both the direction and magnitude of J ("scalar RR"), thereby 
driving stars to near-radial orbits that interact strongly with 
the MBH. RR in a near-spherical potential can only change 
the direction of J ("vector RR"). 

RR is particularly relevant in the potential near a 
MBH, where compact EMRI candidates originate. 

' Present address: Israel Aerospace Industries, MLM Division, PO Box 
45, Beer-Yaakov, 70350, Israel 
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iHopman & Alexander! (|2006|) show that RR dominates 
EMRI source dynamics. Depending on its still poorly- 
determined efficiency, RR can either increase the EMRI rate 
over that predicted assuming NR only, or if too efficient, it 
can strongly suppress the EMRI rate by throwing the compact 
remnants into infall (plunge) orbits (cf Fig. |6]below) that emit 
a single, non-periodic and hard to detect GW burst. A prime 
motivation for the systematic numerical investigation of RR 
efficiency presented here, are the still open questions about 
the implications of RR for EMRI rates and orbital properties. 
This paper is organized as follows. In f|2]we briefly review 
the theory of NR and RR relaxation and derive a new relation 
between scalar and vector RR. In <l3]we describe our method 
of analyzing and quantifying the effects of RR in A^-body sim- 
ulations, which are described in f|4] We present our results in 
^and discuss and summarize them in <l6] 

2. THEORY 

2.1. Non-coherent Relaxation (NR) 

The NR time for iJ-relaxation, Tjfp, corresponds to the 
time it takes non-coherent 2-body interactions to change the 
stellar orbital energy by order of itself, \/S.E\ ^ E (by stel- 
lar dynamical definition convention, E > {) for a bound or- 
bit). Similarly, the NR time for J-relaxation, T-^-^, corre- 
sponds to the time it takes the stellar orbital angular momen- 
tum to change by order of the circular angular momentum 
I A J| - Jc, where near a MBH of mass M, J^ = GM/V2E. 
The _E-relaxation timescale can be estimated by considering 
the rate F of gravitational collisions in a system of size R 
at a relative velocity v, between a test star and N field stars 
of mass m and space density n ^ N/R^, at the minimal 
impact parameter where the small angle deflection assump- 
tion still holds, Tmin ^ Gm/v^. The collision rate is then 
r ^ nvrf^^j^ ^ G'^m^n/v^. Taking into account also col- 
lisions at larger impact parameters increases the rate by the 
Coulomb logarithm factor InA ^ \vi{R/rmin)- Therefore, 
T^^^v^/{G^m^n\nK). Near the MBH w ^ ~ GM/R, and so 
In A '^ In Q, where Q = M/m is the mass ratio. 

When the stars move under the influence of the cen- 



tral MBH (Q 3> N), the relaxation time can be ex- 
pressed as Tjfj^ ~ (M/m)2p/(iVlnA), where P = 
^■n ^ B? I GM is the Keplerian period. Following the nota- 
tion of Rauch & Tremaine ( 1996) (RT96), the NR changes in 
E, J and J over the dimensionless time-lag T={t2 — ti)/Pi 
are 



\AE\ /E= \E2 - Ei\/Ei = aAVN{m/M)V^ , (1) 
|AJ| /Jc = \J2 - Ji\/Jc,i = ?7sA\/iV(m/A'f)V7, (2) 
|AJ| /Jc= |J2 - Jil /Jc.i - r,,^^/N{m|M)^, (3) 

where aA = a^\\\K and rjs.vK = ?/s.t>\/lnA are dimension- 
less constants, whose exact values are system-dependent and 
cannot be estimated with accuracy without detailed calcula- 
tions or simulations. The corresponding NR timescales are 
related to these coefficients by Tjfpj = {M/m)^P/{Na\), 
{M/7nfP/{N7f^^) and T^^ = {M/mfP/{N7i^). 






2.2. Resonant Relaxation (RR) 



When the potential has symmetries that restrict the orbital 
evolution, for example to fixed ellipses in the potential of a 
point mass, or to planar annuli in a spherical potential, the per- 
turbations on a test star are no longer random, but correlated. 
This leads to a coherent changes in J on times P <^t <t^, 
AJ = Tt, by the residual torque |T| -- ^/NGm/R exerted 
by the N randomly oriented, orbit-averaged mass distribu- 
tions of the surrounding stars (mass wires for elliptical orbits 
in a Kepler potential, mass annuli for rosette-like orbits in a 
spherical potential). The coherence time t^ is set by devi- 
ations from the true symmetry, which lead to a gradual or- 
bital drift and to the randomization of T. For example, the 
enclosed stellar mass leads to non-Keplerian retrograde pre- 
cession; General Relativity leads to prograde precession. Ul- 
timately, the coherent torques themselves randomize the or- 
bits (alternatively, this can be viewed as the result of potential 
fluctuations due to the finite number of stars). The effective 
coherence time is set by the shortest de-coherence (quench- 
ing) process in the system. The accumulated change over t^, 
|AJ„| ~ |Ti(j|, then becomes the basic step-size, or mean 
free path in J-space, for the long-term {t ^ t^) non-coherent 
(oc ^/t) relaxation of J. Since this step-size is large, RR can be 
much faster than NR. The RR timescale Trr is then defined 
by AJ/Jc = {AJ^/Jc)^/tJt^= ^t/TRR. Note that the re- 
laxation of E is not affected by RR because the potential of 
the system is stationary on the coherence timescale, and so 
E changes incoherently on all time scales. The torques ex- 
erted by elliptical mass wires in a Kepler potential can change 
both the direction and magnitude of J. In contrast, the torques 
exerted by planar annuli can only change the direction of J. 

Here we consider only Newtonian dynamics. The coher- 
ence timescale for scalar RR is determined by the time it takes 
for the orbital apsis to precess by angle ^ vr due to the poten- 
tial of the enclosed stellar mass ("mass precession"). 



t^=tM = AMiM/Nm)P, 



(4) 



where Am is an 0(1) factor reflecting the approximations 
in this estimate. The coherence timescale for vector RR 
is determined by the time it takes for the coherent torques 
to change A |J| ^ Jc (alternatively, this is the timescale 
tip — (0/^0*)^/2 to accumulate 0(1) fluctuations in the 
stellar potential A0* ^ ^/NGm/R relative to the total gravi- 
tational potential (p as the stars rotate by ^^ tt on their orbits). 



C = ^0 = A^ (VN/2fi) P ~ \A^{M/m) /2Vn] P 



(5) 

where A^f, is an 0(1) factor, iJL — Nm/{M + Nm) and where 
the approximate equality is for the Keplerian limit Nra ^ 
M. Following RT96, the RR changes in J and J during the 
coherent phase {t<t^) can be expressed as 



|AJ| /Jc=\J2 - Ji\ /Jc.i = AViV(m/Af)r , (6) 
|AJ| /J,= IJ2 - Jil /Jc,i = /3„ViV(m/Af)T , (7) 

where the 0(1) dimensionless coefficients jds and (3^ depend 
on the parameters of the system and reflect the uncertainties 
introduced by the various approximations and simplification 
of this analysis. Accurate determination of their values re- 
quires detailed calculations or simulations. 
The scalar RR change A J on time-lags t^tm is then 



\AJ\/J,= \J2-Ji\/JcA=f3s^AM{m/M)V^, (8) 
and the scalar RR timescale is 



T,ij,^[{M/m)/AMP's]P- 



(9) 



The RR efficiency fa c tor y = (/3s//3s,rt96)^ defined by 
iHqpman & Alexandeij (l2006h expresses how much shorter 
-^RR i^ relative to the value estimated by RT96. Scalar RR 
is faster than NR by a factor T^Jr/Irr oc (Af/TO)/iVln A. 
Similarly, the vector RR change | Aj| on time-lags r ^ r^ in 
the Keplerian limit is 



\A3\/J,= \32-Ji\/Jc,i^l3v\I^A^VN{m/M)VT, 



and the vector RR timescale is 



/TlJ 

^RR — 



2iM/m) VNA^P'^ 



(10) 



(11) 



RT96 performed a limited set of near- Keplerian simulations 
to check their predictions. They analyzed the results both star 
by star and in the average and observed the coherent growth of 
AJ/Jc and | AJ| /Jc relative to the simulation's initial values. 
Although the evolution of the of these quantities for any single 
star was very noisy and the proportionality factors had a very 
large scatter, RR was clearly observed, as predicted. 

2.2.1. Relation between scalar and vector RR 

The population averages of the scalar and vector coeffi- 
cients rjs^v and (3s,v, are not independent quantities. Rather, 
{Vv) = c {r]s), and (/3„) — c (/3s) (on timescales t < t^), where 
the constant c depends on the averaging procedure. Here we 
average |AJ| (r) and AJ(r) by the rms over the stellar pop- 
ulation (Sj3]l. We focus on the limit where | AJ| / J <C 1, and 

i2\ 



assume that the change is isotropic on average, ( I AJ| 

3 ( \AJi\ \ = (5^, {i = X, y, z), as is indicated by our simu- 
lations (fQ. The rms of vector RR is 5. Defining the z-axis 
along J, then A J — AJ^ in the small change limit, and its 
rms is (5/\/3- Therefore (see Eqs. |2] |3] IS] 13 ?7t, = \/3?7t, for 



all r and I3y — y/SPs for r <C tm, t^, as is indeed seen in the 
simulations-' (Fig. |2]i. 

The constant ratio between the population averages of the 
vector and scalar coefficients is in part a geometrical effect 
(ID vs 3D changes) and in part a reflection of the symmetries 
of the perturbations (isotropic A J). This strict proportionality 
is not expected to hold when |AJ| / J~ 1, for example when 
only very eccentric stars ( J ^ 0) are included in the sample, 
or for very long time-lags. 

3. RR DETECTION BY AUTO-CORRELATION ANALYSIS 

Given the high computational cost of the A^-body simula- 
tions, and the very large variance in the evolution of individ- 
ual orbits, it is essential to extract the RR signal as efficiently 
and robustly as possible from the simulated data. After some 
experimentation, we adopted the auto-correlation analysis for 
detecting and measuring RR in A^-body simulation snapshots. 
Our procedure is as follows. 

(1) The stellar phase space coordinates are transformed to 
the rest-frame of the MBH, which is almost identical with the 
center of mass in the near-Keplerian system close to the MBH. 

(2) The energy (i^^ ), angular momentum (J-' ), circular 
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angular momentum (J^^ ) and Keplerian period {P^ ) are 
calculated for the n'th star (n = 1 . . . A^) at discrete times ti 
in the simulation ("snapshots"). 
(3) To make full use of the data, we assign a normalized 
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time-lag t - ^ —{tj — ti)/?^ ' for each pair of times {tj >ti). 
For each of the N stars, we calculate the normalized energy 
and angular momentum differences at all lags, (AE/E) 
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(4) The differences from all stars are binned into discrete r- 
bins according to their associated t^j. The bin rms values and 
their standard deviations are calculated and plotted against the 
bin's average time-lag r = {tij), thereby creating the auto- 
correlation curve. 

By using all possible time-lags recorded in the simulation, 
this approach makes maximal use of the entire data set and 
averages over the strongly fluctuating individual relaxation 
curves (cf the RT96 procedure, ii2.2b . However, this method is 
not entirely free of bin-to-bin bias. Since the number of orbital 

periods completed by a star with a mean period P^ over the 
simulation time igim is isim/^^"'' = Tmax, long-period stars 
will not contribute to a high-r bin, (r)^ if Tmlx < (t)^. Con- 
versely, short-period stars will not contribute to a low-r bin. 



if r^L > 



where T^y„ == minj(ij+i -ij)/F(") is 
set by the minimal time difference between consecutive snap- 
shots. In extreme cases some stars may not contribute to the 
relaxation curve at any (r) j. Since long and short period stars 
could well have systematically different responses to RR, this 
introduces bias to the relaxation curve. This bias can be min- 
imized, at the cost of losing some information, by using only 
the middle range of the r -bins, or at a computational cost, by 
longer simulations with a higher snapshot rate. 

4. SIMULATIONS 

Our A^-body code uses a 5th order Runge-Kutta integra- 
tor with individual time-steps and pair-wise K-S regulariza- 

^ It can be shown that if the population's mean absolute difference ( | A J | ) 
is used to define /3s,v, then /3„ =2/3s. 
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phase [rad] 

Fig. 1. — Phase drift in a 2-body system over a time of Tgim = 1-35 X 
10* for a Q = 3 X 10~^ mass ratio and an extreme eccentricity of e = 
0.99995. The line is the initial phase curve (r as function of i/)); the dots 
are the simulated data. The phase drift of Aip = 0.2875 corresponds to 
TA4, c^ 3.3 X 10^. 
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Fig. 2. — The measured (points) and fitted (lines) auto-correlation curves 
for AE/E, I A J| /Jc and | AJ| /Jc in a Q = 10^, 7 = 1.5 simulation with 
A'^ = 200 particles. The points are the results of the A'^-body simulation, the 
thick lines are the predicted theoretical curves; the thin straight lines show 
the asymptotic (t S> tq ) linear behavior. The best fit parameters and related 
quantities are also listed. 



tion ('Kusta anheimo & Stiefell [19651) . without gravity soften- 
ing. The time steps were chosen to conserve total energy at 
the level of A^tot/^tot ^ O{10-^) - (3(10-*), well below 
the NR energy changes expected in the simulations. 

The A^-body code must maintain phase coherence to a high 
enough accuracy to reliably simulate RR. It is thus impor- 
tant to verify that the numerical phase drift in the simulation 
is much smaller than that due to physical processes, such as 
mass precession. Experience shows that most of the phase 
drift Aijj accumulates near periapse, where the acceleration is 
largest, and is therefore more pronounced in eccentric orbits. 
We estimated the phase drift by integrating a highly eccentric 
(e — 0.99995) 2-body system over many orbits, and plotted 
the relative distance between the 2 masses as function of the 
orbital phase (Fig. [Til, as it evolved during the simulation. The 
phase spread relative to the theoretical phase curve, Aijj, was 
measured near ip — tt/2. We define the phase de-coherence 
timescale as ta^ = {Tr/Aip)Tsiai- We estimate conservatively 
that our near-Keplerian simulation typically have ta^ > 10^ 
(This is for highly eccentric orbits. We do not detect any phase 



drift up to T '^ 10^ in orbits with moderate eccentricities). 
Since tav > t^ ^ Tsim for the models simulated here (cf Fig. 
|2ll, it can be safely neglected. 

Our simulations consisted of 200 particles (including the 
MBH as a free particle). The initial orbital semi-major axes 
were randomly drawn from a p{a)da ex a^^'^da distribution 
for 7=1, 1.5, 1.75, for a in the range (0, 1/2), with eccen- 
tricities drawn from a thermal p{e)de = 2ede distribution, 
with random phases, orbital orientations and isotropic veloci- 
ties. This distribution corresponds to an r ' number density 
distribution with an outer cutoff at radius r = R—l from the 
MBH (in dimensionless units where G = 1, M + Nm = 1). 
These stellar cusps span a wide ran ge of possible physical 
scenarios (e.g. iBahcall & WoHl[l977h . and in particular those 
of LISA targets, which are expected to be relaxed galactic nu- 
clei dAlexander 2007). A typical simulation lasted a few x 100 
system orbital times and resulted in few x 100 snapshots of the 
system configuration. In order to decrease the statistical errors 
we ran ngim = 5-8 simulations with different initial conditions 
for each of the (7,Q) models we studied. 

In Fig. (|2|i we plot the measured auto-correlation curves in a 
typical simulation for r<T^: {AE/E){t), (AJ/Jc)(t) and 
(AJ/Jc)(r). The curves reflect the joint effects of NR and 
RR. The coefficients 77^ „ and /3s, i, are measured by fitting the 
auto-correlation curves in the coherent phase to the functions 



I AJ| /J, = ^{m/M)^it^T + PIt 



I A J| /J, = ^{m/M)^ril^T + I31t^ , (12) 

(13) 

where the two terms in the square root express the contri- 
butions of NR and RR, respectively. Note that on short 
timescales, t < tq = {tIsa/Ps)'^, the RR auto-correlation curve 
rises as ^/r due to the effect of NR. It then rises as r in the 
RR-dominated coherent phase at times tq < t < Tuj, before 
turning over again to a y/r rise in the accelerated random- 
walk phase at T > Tj^ (not shown in Fig. |2]i. The NR and RR 
parameters were derived from the best 2-parameterfits of Eqs. 
(fn [T2l[T3] ) to the data. To control the star to star variance, we 
limited our analysis to time-lag bins that sampled at least 0.75 
of the stars in the simulation (fQ. The excellent fit of the data 
points to the predicted auto-correlation curves seen in Fig. (|2|i 
indicates that RR is present and measurable. 

5. RESULTS 

Although the auto-correlation analysis stabilizes against 
star to star scatter in a single simulation, we still find a large 
simulation to simulation scatter in the derived values of the 
coefficients. We therefore constructed a large grid of near- 
Keplerian models (within the computational time limitations), 
where RR should be clearly detected. 

We summarize our results in table ([T]i and in Figs. ([3] ID 
|5]). The coefficients a^ and rj^^yA do not directly express the 

intrinsic properties of NR, since they should vary as VlnA ~ 
a/IiiQ ~ 4(1 ± 0.07) over the Q-range of our models. This 
small fractional difference and the relatively low statistics of 
our simulation suite may explain why we do not detect a clear 
Q-dependence in these quantities. Since we do not see a clear 
7 or Q dependence, we adopt as the best fit estimates of the 
values of a, rjgy and (3gy their grand average over all the 
simulations (table[T]i. The measured ratio {f3v/Ps) = 1-85 is 
consistent with the predicted ratio of \/3— 1-73 ( i)2.2.1| i. 

RT96 derived from their simulations smaller mean values 
for Ps (0.53, as compared to 1.07 here) and for aA and r/sA 
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Fig. 3. — Measured NR energy coefficient qa as function of mass ratio Q 
and for stellar density profiles with logarithmic slopes of 7 = 1, 1.5, 1.75. A 
oc y/lnQ curve is shown to guide the eye. 
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Fig. 4. — Same as Fig. (5), for the measured NR scalar angular mo- 
mentum coefficient r]s\ (bottom points) and the vector angular momentum 
coefficients r]y\ (top points). 
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Fig. 5. — Same as Fig. (5), for the measured RR scalar angular momentum 
coefficient f3s (bottom points) and vector angular momentum coefficient f3y 
(top points). RR does not depend on a Coulomb factor 




efficiency x=(Ps/Ps,rt96)^ 



Fig. 6. — The change in the EMRI and plunge rates relative to that pre- 
dicted assuming only NR, as function of the RR efficiency x (Adapted from 
[Hopman & Alexander 2006, Fig. 5). The RT96 estimate of /Bs (x = 1) pre- 
dicts an increase of X 7. 1 in the EMRI rate due to RR, close to the maximum. 
Our (EKA08) new measured efficient RR (x = 4.1) predicts a higher plunge 
rate and thus a smaller increase of x4.6 in the EMRI rate. 

(3.09 and 1.37 as compared to 14.76 and 4.25 here). At 
least part of the difference in a\ and tjsa can be traced to 
their use of a softening length e = lO^^R ^ Gm/v^ ~ 
{m/M)R in the calculation of the gravitational force. This 
decreases the effective value of the Coul omb factor, since 
A ^ i?/max(e, Gm/w^) = 10^ and so VlnA ~ 2.1, about 
twice as small as in our simulations. Indeed, RT96 noted that 
decreasing the softening length to e = lO^^i? led to an in- 
creased value of QfA = 5.5 ± 0.2. 

It is difficult to trace the specific reason for the discrepancy 
between our best estimate value for Ps and that derived by 
RT96, given the many differences in both the simulations and 
the methods of analysis. Our statistics are better due to the 
larger number of simulations and more efficient use of the 
data, and the analysis here is more rigorous. We therefore 



briefly consi der the imphcations of this revi sed value of /3s for 
EMRI rates. iHopman & Alexanded (l2006l) parametrized the 
RR efficiency by a factor x — [/3s//3s,rt96] . and derived the 
X dependence of the branching ratio of the EMRI and infall 
(plunge) rates (Fig. |6]l. The RT96 value x = 1 happens to lie 
very close to the maximum of the RR-accelerated EMRI rate. 
We find here f3s/Ps,KT96 ^ 2, which corresponds to a factor 
~ 5 increase in the EMRI rate compared to that estimated 
for NR only, but is a factor ^1.5 smaller than implied by 
the RT96 value, because the higher RR efficiency leads to a 
higher plunge rate at the expense of the inspiral rate. 

6. DISCUSSION AND SUMMARY 

We characterized and measured the mean efficiency coeffi- 
cients of NR (ofA, ?7s,A, '7)j,a) and RR {(3s, (3y) in Newtonian 
A^-body simulations of isotropic, thermal, near-Keplerian 
stellar cusps around a MBH. We derived a simple analytical 
form for the rms RR auto-correlation curves of J and J, and 
showed that the two are simply proportional to each other. We 
then measured these coefficients in a large suite of small scale 
A^-body simulations with different stellar density distributions 
and MBH/star mass ratios. We found no statistically signif- 
icant trends in the values of these coefficients as function of 
the system properties. This may require better statistics. 

Our measured RR efficiency suggests that RR increases 
the EMRI rate by a factor of ^ 5 above what is predicted 
for NR only. This estimate of RR efficiency is consistent 
with that suggested by the analysis of the dynamical prop- 
erti es of the different stellar popu lations in the Galactic Cen- 
ter (iHopman & Alexandedl2006h . However, this conclusion 
is still preliminary, since several important open issues re- 
main, which should be addressed by larger scale simula- 
tions. These include (1) The dependence of the RR coeffi- 
cie nts on the orbit of the test star, for example its eccentric- 
ity dGiirkan & HopmanlllOOTh . We find that A^ - 200 is not 
enough for reliable statistics on sub-samples within given E 
or y-bins. The eccentricity dependence of RR is particularly 
relevant for the supply rate of stars to the MBH from J ^ or- 
bits (the loss-cone refilling problem). (2) The effects of a stel- 
lar mass spectrum. This will likely affect RR by decreasing 
the RR timescale, and changing the stellar density distribution 
through strong mass segregation (Alexander 2007; Alexander 
& Hopman 2008, in prep.). (3) The robustness of RR against 
perturbations from the larger non-Keplerian stellar system in 
which the inner near-Keplerian region of interest is embed- 
ded. (4) The role of post-Newtonian effects in RR, such as 
General Relativistic precession and GW emission. These are 
expected to play a key role in enabling inspiral by quench- 
ing RR just as the compact remnant e nters the EMRI phase, 
and in regulating the GW inspiral rate dHopman & Alexandeij 
l2006h . 
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